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ABSTRACT 

We observed Sgr A* and two extragalactic radio sources nearby in angle with the VLBA 
over a period of two years and measured relative positions with an accuracy approaching 0.1 
mas. The apparent proper motion of Sgr A* relative to J1745-283 is 5.90 ± 0.4 mas y~^, almost 
entirely in the plane of the Galaxy. The effects of the orbit of the Sun around the Galactic 
Center can account for this motion, and any residual proper motion of Sgr A*, with respect to 
extragalactic sources, is less than about 20 km s~^ . Assuming that Sgr A* is at rest at the 
center of the Galaxy, we estimate that the circular rotation speed in the Galaxy at the position 
of the Sun, Gq , is 219 ± 20 km s^^ , scaled by i?o/8.0 kpc. 

Current observations are consistent with Sgr A* containing all of the nearly 2.6 x 10^ Mq, 
deduced from stellar proper motions, in the form of a massive black hole. While the low 
luminosity of Sgr A*, for example, might possibly have come from a contact binary containing 
of order 10 M© , the lack of substantial motion rules out a "stellar" origin for Sgr A*. The very 
slow speed of Sgr A* yields a lower limit to the mass of Sgr A* of about 1,000 M© . Even for 
this mass, Sgr A* appears to be radiating at less than 0.1% of its Eddington limit. 

Subject headings: Individual Sources: Sgr A*; Black Holes; Galaxy: Center, Fundamental 
Parameters, Structure; Astrometry 

1. Introduction 

Sgr A* is a compact radio source, similar to weak active nuclei found in other galaxies. Since 
its discovery more than two decades ago by Balick & Brown (1974), the possibility that Sgr A* is a 
super-massive (~ 10^ M©) black hole has been actively considered. However, its radio luminosity of 
w 10^ L© (Serabyn et al. 1997) and its estimated total luminosity of 10^ L© are many orders of magnitude 
below that possible from a ~ 10^ M© black hole. Thus, on the basis of its spectral energy distribution 
Sgr A* could be an unusual contact binary containing ^ 10 M© and radiating near its Eddington limit. 

Recently, Eckart & Genzel (1997) and Chez et al. (1999) measured proper motions of stars near the 
position of Sgr A*, as determined by Menten et al. (1997). Stellar speeds in excess of 1000 km s~^ at a 
distance of 0.015 pc from Sgr A* indicate a central mass of 2.6 x 10^ M©. While these dramatic results 
are consistent with the theory that Sgr A* is a super-massive black hole, it is still conceivable that most of 
the central mass could come from a combination of stars and perhaps some form of dark matter. Clearly, 
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independent constraints on the mass of Sgr A* are needed to establish whether or not it is a super-massive 
black hole nearly at rest at the dynamical center of the Galaxy. 

The apparent motion of Sgr A* can be used to estimate the mass and elucidate the nature of this 
unusual source. An apparent motion of Sgr A* can be attributed to at least three possible components: 1) 
a secular motion induced by the orbit of the Sun about the Galactic Center, 2) a yearly oscillation owing to 
the Earth's orbital motion around the Sun (trigonometric parallax) , and 3) a possible motion of Sgr A* with 
respect to the dynamical center of the Galaxy. Measurement of, or limits for, these components of Sgr A*'s 
apparent motion can provide unique information on the circular rotation speed {Qq ) of the Local Standard 
of Rest (LSR) and the peculiar motion of the Sun {Vq ) , the distance to the center of the Galaxy {Rq ) , and 
the nature of Sgr A* itself. 

In 1991 we started a program with the Very Long Baseline Array (VLBA) of the National Radio 
Astronomy Observatory |^ (NRAO) to measure the apparent motion of Sgr A* . In principle, VLBA 
observations of Sgr A* , phase-referenced to extragalactic radio sources, can achieve an accuracy sufhcient 
to detect secular motions of ^ 1 kms^^ for a source at the distance of the Galactic Center. However, 
achieving this accuracy is quite challenging technically as it involves observing at short wavelengths (7 mm) , 
in order to minimize the effects of interstellar scattering toward Sgr A*, phase-referencing to extragalactic 
sources, and careful modeling of atmospheric effects (because of the low source elevations). 

In the early years of the project, we searched for strong, compact, extragalactic sources nearby in 
angle to Sgr A*and worked toward an optimum observing strategy. In this paper, we report results of 
observations spanning a two year period from 1995 to 1997. Our VLBA images clearly show movement of 
Sgr A* with respect to extragalactic sources over many synthesized beams. While the current positional 
accuracy is inadequate to determine a trigonometric parallax, the secular motion of Sgr A* is easily 
measured. This yields an accurate estimate of the angular rotation rate of the Sun around the Galactic 
Center, (8o + Vo)/i?o, and places interesting limits on the mass of the black hole candidate responsible for 
the radio emission. Our results, and those of Backer and Sramek (1999) from observations with the Very 
Large Array (VLA) over a 15 year period, indicate that the apparent proper motion of Sgr A* is dominated 
by the orbit of the Sun about the Galactic Center and that any peculiar motion of Sgr A* is very small. 



2. Observations 

Our successful observations using the VLBA were conducted in the late-night, early- morning periods 
of 1995 March 4, 1996 March 20 and 31, and 1997 March 16 and 27. (Observations attempted during two 
August evenings in 1996 experienced high water vapor turbulence in the atmosphere, and phase referencing 
was not successful.) The observing frequency was 43.2 GHz and we observed four 8 MHz bands, each at 
right and left circular polarization. We employed 2-bit sampling at the Nyquist rate, which required the 
maximum aggregate sampling rate supported by the VLBA of 256 megabits per second. Only the inner five 
VLBA stations (Fort Davis, Los Alamos, Pie Town, Kitt Peak, and Owens Valley) were used, as baselines 
longer than about 1500 km heavily resolve the scatter broadened image of Sgr A* at 43 GHz (eg. Bower & 
Backer 1998). 

The observing sequence involved rapid switching between compact extragalactic sources and Sgr A*. 



^The National Radio Astronomy Observatory is operated by Associated Universities Inc., under cooperative agreement with 
the National Science Foundation. 
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Two sources, J1745-283 and J1748 291, from the catalog of Zoonematkermani et al. (1990) were found to 
be strong enough (> 10 mJy at 43 GHz) to serve as reference background sources. These sources are two of 
the three used by Backer & Sramek (1999) in their program to measure the proper motion of Sgr A* with 
the VLA; their third background source, J1740-294, proved to have a steep spectrum and was too weak for 
inclusion in our 43 GHz program. Wc switched among the sources repeating the following pattern: Sgr A*, 
J1745-283, Sgr A*, J1748-291. Sources were changed every 15 seconds, typically achieving 7 seconds 
of on-source data, except for the earliest observation in 1995 when we were experimenting with longer 
switching times. We used Sgr A* as the phase-reference source, because it is considerably stronger than 
the background sources and could be detected on individual baselines with signal-to-noise ratios typically 
between 10 and 20 in the 7 seconds of available on-source time. 

We edited and calibrated data using standard tasks in the Astronomical Image Processing System 
(AIPS) designed for VLBA data. This involved applying data flagging tables generated by the on-line 
antenna and correlator systems, station gain curves, and system temperature measurements. We solved 
for station-dependent, intermediate- frequency band delays and phases on a strong, compact source 
(NRAO 530). After applying these corrections, the multi-band data for Sgr A* could be combined 
coherently and interferometer phases as a function of time determined. These phase solutions were 
examined by an AIPS task specially written for our observations that looked for and flagged data when 
baseline-dependent phases on adjacent Sgr A* scans changed by more than one radian. Under good 
weather conditions between 10 and 30% of the data were discarded by this process. This provided relatively 
unambiguous "phase connection" for the remaining data and allowed removal of most of the effects of 
short-term atmospheric fluctuations from all sources. (We note that during average-to-poor weather 
conditions, our phase measurements on Sgr A* every 30 seconds were not frequent enough to provide 
unambiguous phase connection. Thus, our 15 second switching time is probably an optimum trade-off 
between on-source duty cycle and atmospheric coherence losses.) Data calibrated in this manner produced 
high (eg, 50:1) dynamic range maps of all sources with little or no spatial blurring. The images of the 
background sources appeared less resolved than that of Sgr A*, with no signs of complex or multiple 
component structures. 

We found that the differences in relative positions between a background source and Sgr A* for closely 
spaced (w 10 d) epochs were ~ 1 mas. These differences exceeded the formal precision, estimated by the 
least-squares fitting process, typically by a large factor. Since the observational conditions and data analysis 
were nearly identical for these epochs, small geometric errors (eg, in baselines, source coordinates, or Earth's 
orientation parameters) are unlikely to yield position shifts of this magnitude given the small angular 
separation of Sgr A* and the background sources. Therefore, we evaluated the possibility that refractive 
scattering of the radio waves in the interstellar medium or modeling errors for the Earth's atmospheric 
propagation delay could be responsible. 

Refractive scattering can cause changes in the apparent flux density and position of a radio source. 
Gwinn et al. (1988) published VLBI observations that limit refractive position wander for water vapor 
masers in Sgr B2N, a star forming region close to the Galactic Center. For these maser spots, their 22 GHz 

observations revealed a diffractive scattering size of 0.3 mas and an upper limit to a Gaussian component of 
refractive position wander of 0.018 mas. Theoretical estimates of refractive position wander, based on the 
diffractive scattering size, by Romani, Narayan & Blandford (1986) agree with this limit for a Kolmogorov 
electron density spectrum. Assuming that the refractive effects scale as the diffractive scattering size, we 
expect any refractive wander of Sgr A* at 43 GHz to be < 0.04 mas, a value about a factor of 25 smaller 
than our observed position differences. 
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After careful study of the data, we concluded that the most likely source of relative position error is a 
small error in the atmospheric model used by the VLBA correlator. The following simple analysis supports 
this view: The phase-delay of the neutral atmosphere, r, can be approximated by Tq sec Z, where Tq is 
the vertical phase-delay and Z is the local source zenith angle. When measuring the difference in position 
of two sources separated in zenith angle by AZ, a first-order Taylor expansion of r yields the expected 
differenced phase-delay error for a single antenna: 

dT 

At w = rosecZtanZAZ . (1) 

oZ 

The seasonally-averaged atmospheric model (Niell 1996) used by the VLBA correlator is likely to 
miss-estimate tq by about 0.1 nsec, equivalent to a zenith phase-delay of ~ 3 cm in path length. This comes 

mostly from the highly variable contribution by water vapor (eg, Treuhaft & Lanyi 1987). Based on Eq. (1), 
this should result in an antenna-dependent error of At ~ 0.3 cm for typical source zenith angles of « 70° 
and for our typical source separations of AZ w 0.012 rad(R:! 0.7°). Since atmospheric errors are largely 
unc:orrelated for different antennas, on an interfc;romctcr baseline at an observing wavelength of 0.7 cm, we 
would expect a relative position shift of roughly \/2 x 0.3 cm or w 70% of a fringe spacing. This corresponds 
to ~ 0.4 and 1.7 mas in the easterly and northerly directions, respectively, for our longer baselines. This 
effect should only partially cancel among the different baselines and can explain the position errors seen in 
the raw maps made from observations made w 10 days apart. 

In order to improve our relative position measurements, we modeled simultaneously our differenced- 
phase data for the "J1745 283 minus Sgr A*" and "J1748 291 minus Sgr A*" source pairs. The model 
allowed for a relative position shift for each source pair and a single vertical atmospheric delay error in 
the correlator model for each antenna. This approach significantly improved the accuracy of the relative 
position measurements as evidenced by the smaller deviations in relative positions for observations closely 
spaced in time. The vertical atmospheric delay parameters typically indicated a correlator model error of 
a few cm and these parameters were estimated with uncertainties of about 1 cm. Using this approach, we 
would estimate from Eq. (1) and the above discussion that relative position errors should be w 0.1 and 
pa 0.4 mas in the East-West and North-South coordinates, respectively, for one day's observation. 

The differenced-phases often displayed post-fit residuals of ~ 30 degrees of phase, which were correlated 
over periods of hours. Assuming equal and uncorrclatcd contributions from the two antennas forming an 
interferometer pair, this suggests a delay change of about 0.1 nsec of time or about 3 cm of uncompensated 
path length. Since typical source zenith angles were about 70 degrees, this corresponds to a vertical delay 
change in the atmosphere above each antenna of about 0.03 nsec or about 1 cm of path length. This 
behavior is consistent with expected large-scale changes in the atmospheric delay, and it suggests that 
significant improvement can be obtained by monitoring and correcting for large scale atmospheric changes. 

The data in Table 1 summarize our relative position measurements. The data taken on 1995 March 4 
were of poor quality, only the stronger of the two background sources (J1745-283) was detected, and the 
positional accuracy was significantly worse than for the later epochs. Positions of the strongest background 
source, J 1745-283, phase-referenced to Sgr A*, for epochs spanning 2 years are plotted in Fig. 1 with open 
circles in the sense Sgr A* relative to J1745-283. They indicate a clear apparent motion for Sgr A* relative 
to J 1745-283, consistent in magnitude and direction with the reflex motion of the Sun around the Galactic 
Center (sec §3.1). The positions in the East- West direction have typical uncertainties of about 0.1 mas, as 
estimated from the scatter of the post-fit position residuals about a straight-line motion. It is interesting 
to note that, while it takes w 220 My for the Sun to complete an orbit around the Galactic Center, the 
East-West component of the parallax from only 10 days motion can be detected with the VLBA! The 
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position uncertainties in the North-South direction are larger, about 0.4 mas, owing to the low declination 
of the Galactic Center. 

The apparent motions of Sgr A* relative to J1745-283 over a 2 year time period and .J1748 291 
relative to J 1745-283 over a 1 year time period are given in Table 2. The uncertainties in Table 2 include 
estimates of the systematic effects, dominated by errors in modeling of atmospheric effects, as discussed 
above. Assuming that J1745-283 is sufficiently distant that it has negligible intrinsic angular motion, 
Sgr A*'s apparent motion is — 3.33± 0.1 and — 4.94± 0.4 mas y~^ in the easterly and northerly directions, 
respectively. This motion is shown by the solid lines in Fig. 1. 

As a check on the accuracy of our measurements, we measured relative positions between the 
two calibration sources. These positions are plotted in Fig. 1 (crosses) in the sense "J1748-291 minus 
J1745-283," offset to fit the plotting scale for the "Sgr A* minus J1745-283" data. The best fit motions are 
+0.17 ± 0.14 and — 0.22 ± 0.56 mas y"'^ in the easterly and northerly directions, respectively, as indicated 
with the dashed lines in Fig. 1. The uncertainty in the relative motion of J1748-291 with respect to 
J1745-283 is Ri40% larger than for the motion of Sgr A* with respect to J1745-283, because the angular 
separation of the two background sources (« 1.0 deg) is greater than between Sgr A* and either of the 
background sources (« 0.7 deg). Thus, the background sources display no statistically significant motion 
relative to each other, as expected for extragalactic sources. 

Finally, we have determined the position of Sgr A* relative to an extragalactic source with high 
accuracy and, therefore, can derive an improved absolute position of Sgr A*. VLBI observations carried out 
by the joint NASA/USNO/NRAO geodetic/astrometric array (Eubanks, private communication) detected 
J1745-283 at 8.4 GHz and determined its position in the U.S. Navy 1997-1998 reference frame to be 

J1745-283 a(2000)=17 45 52.4968, (5(2000)=-28 20 26.294 , 
with a uncertainty of about 12 mas. Assuming this result, 
we find the position of Sgr A* measured at 1996.25 to be 

Sgr A* a(2000)=17 45 40.0409, (5(2000)=-29 00 28.118 . 

The uncertainty in this position is dominated by that of J 1745-283. Note that were one not to correct for 
the "large" apparent proper motion of Sgr A*, the position of Sgr A* determined for observations made 
more than 2 years from 1996.25 would be shifted by an amount greater than the « 12 mas uncertainty. 

3. Discussion 

The apparent motion of Sgr A* with respect to background radio sources can be used to estimate the 
rotation of the Galaxy and any peculiar motion of the super-massive black hole candidate Sgr A*. In Fig. 2 

we plot the change in apparent position on the plane of the sky of Sgr A* relative to J 1745-283. The dotted 
line is the variance-weighted least-squares fit to the data, and the solid line denotes the orientation of the 
Galactic Plane. Clearly the apparent motion of Sgr A* is almost entirely in the Galactic Plane. Thus, it 
is natural to convert the apparent motion from equatorial to galactic coordinates. For Sgr A* relative to 
J1745-283, this yields an apparent motion of —5.90 ± 0.35 and +0.20 ± 0.30 mas y^^ in galactic longitude 
and latitude, respectively. The apparent motion in the plane of the Galaxy should be dominated by the 
effects of the orbit of the Sun around the Galactic Center, while the motion out of the plane should contain 
only small terms from the Z-component of the Solar Motion and a possible motion of Sgr A*. In the 
following subsections, we investigate the various components of the apparent motion of Sgr A*, place limits 
on any offset of Sgr A* from the dynamical center of the Galaxy, derive limits on the mass of Sgr A*, and 
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constrain the distribution of dark matter in the Galactic Center. 



3.1. Motion of Sgr A* in the Plane of the Galaxy 

Assuming a distance of 8.0 ± 0.5 kpc (Reid 1993), the apparent angular motion of Sgr A* in the plane 
of the Galaxy translates to —223 ± 19 km s""*^ . The uncertainty includes the effects of measurement errors 
and the 0.5 kpc uncertainty in Rq . Provided that the peculiar motion of Sgr A* is small (see §3.2), this 
corresponds to the reflex of true orbital motion of the Sun around the Galactic Center. This reflex motion 
can be parameterized as a combination of a circular orbit (i.e., of the LSR) and the deviation of the Sun 
from that circular orbit (the Solar Motion). The Solar Motion, determined from Hipparcos data by Dehnen 
& Binney (1998), is 5.25 ±0.62 km s""'^ in the direction of galactic rotation. Removing this component of the 
Solar Motion from the reflex of the apparent motion of Sgr A* yields an estimate for Go of 218 ± 19 km s^^ 
. This value is consistent with most recent estimates of about 220 km s^^ (Kerr & Lynden-Bell 1986) and 
can be scaled for different values of the distance to the Galactic Center by multiplying by i?o/8 kpc. 

The most straightforward comparison of our direct measurement of the angular rotation rate of the 
LSR at the Sun (Qq / Rq ) can be made with Hipparcos measurements based on motions of Ccphcids. 
Feast & Whitelock (1997) conclude that the angular velocity of circular rotation at the Sun, Oq /Rq (= 
Oort's A-B), is 27.19 ± 0.87 km s^^ kpc^i (218 ± 7 km s^^ for Rq = 8.0 kpc). Our value of Bq/Ro, 
obtained by removing the Solar Motion in longitude from the refl,ex of the motion of Sgr A* in longitude, 
is 27.2 ± 1.7 km s"""^ kpc"-'^. The VLBA and Hipparcos measurements are consistent within their joint 
errors, and both measurements are insensitive to the value of Rq , as it is only used to remove the small 
contribution of the Solar Motion. It is important to note that our value of Go /i?o is a true "global" measure 
of the angular rotation rate of the Galaxy. The consistency of the local (A-B) and global measures of 
Qo/Ro suggests that local variations in Galactic dynamics (dBo/dRo) are less than the joint uncertainties 
of about 2 km s^^ kpc"-'^. 

After removing the best estimate of the motion of the Sun around the Galactic Center, our VLBA 
observations yield an estimate of the pecuhar motion of Sgr A* of 0.0 ± \/0.87^ + 1.7^ km s~^ kpc~^ or 
± 15 km s~^ towards positive galactic longitude. This estimate of the "in plane" motion of Sgr A* comes 
from differencing two angular motions. Since this difference is negligible, the uncertainty in Rq does not 
affect this component of the peculiar motion of Sgr A*. Given the excellent agreement in the global and 
local measures of the angular rotation rate of the Galaxy, and the lack of a detected peculiar motion for 
Sgr A*, it is likely that Sgr A* is at the dynamical center of the Galaxy. 

3.2. Motion of Sgr A* out of the Plane of the Galaxy 

Whereas the orbital motion of the Sun (around the Galactic Center) complicates estimates of the "in 
plane" component of the peculiar motion of Sgr A*, motions out of the plane are simpler to interpret. One 
needs only to subtract the small Z-component of the Solar Motion from the observed motion of Sgr A* to 
estimate the out-of-plane component of the peculiar motion of Sgr A*. An implicit assumption in this 
procedure is that the Solar Motion reflects the true peculiar motion of the Sun. Since most estimates of 
the Solar Motion are relative to stars in the solar neighborhood, this assumes that "local" and "global" 
estimates of the Solar Motion are similar. This procedure could be compromised slightly were the solar 
neighborhood to have a significant motion out of the plane of the Galaxy, owing, for example, to a galactic 
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bending or corrugation mode. 

One way to limit the magnitude of a possible difference between a local and global estimate for the 
Solar Motion is to compare motions based on nearby stars with those based on much more distant stars. 
Using stars within about 0.1 kpc Dehnen & Binney (1998) find the Z-component for the Solar Motion to be 
7.17 ± 0.38 kms"\ while Feast & Whitelock (1997) determine a value of 7.61 ± 0.64 kms"^ for stars with 
distances out to a few kpc. Since these values agree within their joint uncertainties of about 0.74 km , it 
seems unlikely that local and globlal values for the Solar Motion could differ by more than about 1 km s~^ . 

The Hipparcos measurements of large numbers of stars in the solar neighborhood provide an excellent 
reference for determining the local solar motion. We adopt the value of Dehnen & Binney (1988), which 
comes from the velocities of more than 10,000 stars within about 0.1 kpc of the Sun. Removing 7.17 km 
s"""^ from our measured apparent motion of Sgr A* out of the plane of the Galaxy, we estimate the peculiar 
motion of Sgr A* to be 15 ± 11 km s^^ toward the north galactic pole (see Table 3). The uncertainty is 
dominated by our proper motion measurements and can be greatly improved by future measurements. Note, 
for example, that increasing the weight (decreasing the estimated uncertainty) of the 1995 measurement 
would decrease the magnitude of our peculiar motion estimate. We do not consider our estimate of the 
peculiar motion of Sgr A* out of the plane of the Galaxy to be statistically distinguishable from a null 
result. 



3.3. Limits on the Mass of Sgr A* 

Our estimates of a peculiar motion of Sgr A* provide an upper limit of about 20 km s~^ each for 
motions in and out of the galactic plane. Since stars in the inner-most regions of the central cluster move 
at speeds in excess of 1000 km s~^ (Eckart & Genzel 1997, Ghez et al. 1999), a central "dark mass" of 
approximately 2.6 x 10^ Mq contained within 0.015 pc of Sgr A* seems required. It is likely, but unproven, 
that most of this mass is contained in a super- massive black hole: Sgr A*. Given the fact that independent 
measurements (Backer & Sramek 1999, and this paper) show that Sgr A* moves at least two orders of 
magnitude slower than its surrounding stars, Sgr A* must be much more massive than the ^ 10 Mq stars 
observed in the central cluster. (See also Gould & Ramirez [1998] for discussion of the implications of a 
lack of apparent acceleration of Sgr A*.) In this section we derive a lower limit to the mass of Sgr A* and 
constrain possible distributions of dark matter, not in the form of a super-massive black hole. 

3.3.1. Virial Theorem 

Unfortunately, the Virial theorem is of little help in relating the masses and velocities of stars to that 
of a central massive black hole. For Virial equilibrium, 

Ts+nh = -^{Us + Ubh) , (2) 

where T and U correspond to the kinetic and potential energy terms and the subscripts s and bh identify 
those associated with the stars and a central, massive, black hole, respectively. Essentially all the kinetic 
energy can be tied up in the stars (Tg) and all the gravitational potential energy found associated with 
the black hole (Ubh)- In this case, attempts to estimate the kinetic energy of the black hole (Tb/j) require 
differencing two large and uncertain quantities and will be essentially useless. 
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3.3.2. Equipartition of Kinetic Energy 

Upper limits on the motion of Sgr A* have been used to infer lower limits on the mass of Sgr A* by 
assuming equipartition of kinetic energy (eg, Backer 1996, Genzel et al. 1997). This is reasonable for stellar 
systems such as globular clusters, where massive stars are found concentrated toward the cluster center and 
move more slowly than lower mass stars. Similar results might also hold for a system involving a central 
black hole and a surrounding stellar cluster, provided the core mass of the cluster greatly exceeds that of the 
black hole. However, given the likely mass dominance of Sgr A* over the stars within 0.015 pc, where high 
stellar speeds have been measured, equipartition of kinetic energy may be an unreliable approximation. 

Indeed, our solar system may prove a better "scale model" (with planets corresponding to stars and 
the Sun corresponding to Sgr A*). The Sun orbits the barycenter of the solar system, approximately in 
a binary orbit with Jupiter. Neglecting small perturbations from other planets, for a binary orbit in the 
center of mass frame, momentum conservation requires that mjvj = Mi^^y^;, where the subscripts J and 
refer to Jupiter and the Sun, respectively. The ratio of the kinetic energy of Jupiter to the Sun is given by 
ttljv'^j/MqVq = MqIvhj. Hence, the kinetic energy of Jupiter exceeds that of the Sun by a factor equal to 
the inverse of the ratio of their masses and equipartition of kinetic energy does not apply. 



3.3.3. Case 1: MsgrA* ~ 2.6 x 10^ Mq 

In order to better evaluate how an upper limit to the motion of Sgr A* can be used to provide a lower 
limit to its mass, we carried out N-body simulations of stars orbiting about a massive black hole. We used 
a simple, direct integration code (NBODYO) of Aarseth (1985), documented by Binney & Tremaine (1987) 
and modified for our purposes. Initial simulations used 255 stars orbiting a 2.6 x 10^ Mq black hole. The 
stellar masses were chosen randomly to represent the upper end of a stellar mass function with a power 
law distribution from 20 down to 2 Mq . The number of stars and their masses used in this simulation are 
comparable to those observed within the central 0.5 arcsec, or 4000 AU (Genzel et al. 1997). Stellar orbits 
were chosen by randomly assigning a distance from the black hole, uniformly distributed in the range 10 to 
10,000 AU, calculating a circular orbital speed, and then adjusting the speed randomly by between ± 20% 
of the circular speed for each of the three Cartesian coordinates. Before starting the N-body integrations, 
the orbital orientations were randomized by rotating the coordinates (and velocity components) through 
three Euler rotations with angles chosen at random. 

The N-body simulations show quasi-random motions of the massive black hole. After relatively short 
periods of time (<C 10, 000 years) a "steady state" condition appeared to be reached. The speed of a typical 
star was about 700 km s~^ at an average distance of 6,000 AU. The motion of Sgr A* changed completely 
in all three coordinates on time scales <C 100 years and was typically ^ 0.1 km s~^ in each coordinate. 
The rapid, but bounded, changes in the motion of Sgr A* suggests that close encoimters with individual 
stars are responsible for most of the observed motions. Assuming this is the case, one can make a simple 
analytical estimate of the expected motion of Sgr A*, owing to close encounters with stars in the dense 
central cluster. 

For a two-body interaction conserving momentum and viewed in the center of mass frame, 

mv = MV , (3) 

where m and v are the mass and speed of the star and M and V are the mass and speed of the black hole, 
respectively. For the case of interest where M ^ m, the orbital speed of the star at periastron, Vp, is given 
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by the well known relation 




GM /1 + eN 



a \l-e) 



(4) 



where G is the gravitation constant, a is the stellar semi-major axis, and e is the orbital eccentricity. 
Defining Vp as the speed of the black hole at periastron, combining Eqs. (3) and (4) yields 



Our observations are only sensitive to orbital periods longer than of order 1 year. Such orbital periods 
occur for stellar semimajor axes greater than about 50 AU (1000 Schwarzschild radii) for a 2.6 x 1O®M0 
black hole. Thus, for our application reasonable values for the parameters in Eq. (5) are as follows: 
m 10 Mq, e ~ 0.5, and a ~ 50 AU, and the expected orbital speed of Sgr A* would be « 0.03 km s~^ . 
(We adopt the periastron speed, instead of the lower average orbital speed, because the influence of many 
orbiting stars will likely increase the speed of Sgr A*, compared to the single star result.) This speed is well 
below our current limit for the motion of Sgr A*. Thus, the simplest interpretation consistent with the fast 
stellar motions and the slow Sgr A* motion is that Sgr A* is a super-massive black hole. 



While the simplest interpretation is that Sgr A* is a 2.6 x 10® Mq black hole, our upper limit on any 
peculiar motion for Sgr A* currently is two to three orders of magnitude above its expected motion for that 
mass. Thus, it seems reasonable to investigate the possibility that the mass within a radius of 0.01 pc is 
not dominated by Sgr A*, but instead is in some form of "dark" matter. In this case, Sgr A* will react 
to the gravitational potential and orbit the center of mass of this dark matter. Below we show that the 
upper limits on the motion of Sgr A* are complimentary to the stellar proper motion results and strongly 
constrain both the mass of Sgr A* and any possible configuration of matter within the central 0.01 pc. 

Fig. 3 displays the enclosed mass versus radius for four mass models that are consistent with the stellar 
proper motions (cf., Genzel et al. 1997, Ghez et al. 1999). These models all yield flat "enclosed- mass versus 
radius" relations at distances ^ 0.01 pc from the center of mass of the system where measurements exist. 
The most centrally condensed mass distribution, a point mass, is shown as the horizontal dash-dot line 
labeled "a" . The least centrally condensed mass distribution plotted is for a Plummer density distribution, 
where density, p, is given by p = po (l + (r/ro)^) for po = 6 x 10"'^'^ Mq pc~'^, tq = 0.01 pc, and a = 5. 
This distribution is shown with a curved dash-dot line labeled "d" . It is difficult to make a physically 
reasonable mass distribution that is significantly less centrally condensed than this and consistent with 
the stellar motion data. These two "extreme" model distributions approximately bound all allowed mass 
distributions; two particular examples of intermediate models are shown in Fig. 3 with dashed curves 
labeled "b" and "c". 

Assuming that Sgr A* has a mass <C 10® M© , it will orbit about the center of mass of the system. 
Orbital speed for a body at a radius, r, from the center of mass is given hy V = y/GMgnci/r, where Mend 
is the enclosed mass at that radius. Setting V = 20 km s~^ , our limit for the motion of Sgr A*, produces 
the sloping solid line in Fig. 3. Only enclosed masses below that line are permitted by our observations. For 
radii greater than about 3 x 10~^ pc (6 AU), this limit rules out all but the least centrally condensed mass 
models. For radii less than about 3 x 10~^ pc, the orbital motion of Sgr A* produces angular excursions 
less than 0.8 mas. In this case, while the orbital speed might greatly exceed 20 km s~^ , we may have failed 




(5) 



3.3.4. Case 2: MggrA* < 2.6 x 10® M, 



■© 
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to detect these excursions owing to our poor temporal sampling and ss 0.4 mas errors in the North-South 
direction. Thus, we do not extend the motion limit line below 3 x 10~^ pc, and at this radius we replace 
the motion limit with a vertical line in Fig. 3. 

The stellar proper motions and the limits on the proper motion of Sgr A* combine to exclude almost 
all of "parameter space" for models of the density distribution of material in the inner 0.1 pc of the Galactic 
Center. The stellar motions exclude "soft" gravitational potentials (i.e., the least centrally condensed mass 
distributions) and the motion limit for Sgr A* excludes "hard" gravitational potentials. Continued VLBA 
observations of Sgr A* over the next five years could reduce the uncertainty in the peculiar motion of 
Sgr A* to about 2 km s~^ (dotted line in Fig. 3) out of the plane of the Galaxy. Improved accuracy for 
positions in individiial VLBA observations, necessary for a trigonometric parallax of Sgr A*, could move 
the small angular excursion limit to < 0.2 mas. This would further and drastically restrict the range of 
possible models for a dominant central dark matter condensation. 

The point-like mass distribution labeled "a" in Fig. 3 essentially requires a super-massive black hole. 
Since we have assumed in this section that this is not Sgr A*, we would be left with the question of why 
a low-mass Sgr A* radiates far more than a super-massive black hole in essentially the same environment. 
Other mass models, and especially those with the most centrally condensed mass distributions (eg, labeled 
"b" in Fig. 3) require exceedingly high mass densities. For example, model "b" has Ri 5 x 10^ Mq within a 
radius of 6 AU, resulting in a density of 10^^ Mq pc~^. Theoretical arguments suggest that such models 
are unlikely to be stable for even 10'^ y, regardless of the composition of the matter (Maoz 1998). 

The hardest cases to exclude on either observational or theoretical grounds are given by the mass 
distributions that are the least centrally condensed, but still consistent with the stellar motion data. For 
this case, one can ask the question, at what mass will perturbations by stars in the central cluster lead to 
detectable motion of Sgr A*? In order to answer this question, we modified the N-body code described 
in § 3.3.3 to include a fixed gravitational potential appropriate for a Plummcr law mass distribution with 
a = 5, po = 6 X 10^^ Mq pc~^, and ro = 0.01 pc (model "d"). Sgr A* was assigned a mass (< 10^ M© ) and 
the entire system, including 254 stars was allowed to evolve in time. For the softer allowable gravitational 
potentials (eg, "c" and "d") there is little enclosed mass within ^ 10^^ pc to bind Sgr A*. We found that 
when Sgr A*'s mass was less than w 3, 000 Mq , Sgr A* gradually moved outward from the center of the 
gravitational potential and achieved orbital speeds in excess of 20 km s~^ . Thus, we conclude that the lack 
of detectable motion for Sgr A* places a conservative lower limit of about 1,000 Mq for the mass associated 
with Sgr A*. 

4. Conclusions 

The apparent proper motion of Sgr A*, relative to extragalactic sources, is consistent with that expected 
from the Sun orbiting the center of the Galaxy. Thus, Sgr A* must be very close to, and most likely at, the 
dynamical center of the Galaxy. In this case, the proper motion measurement gives the angular rotation 
speed at the Sun, (Bq + V0)/i?o, directly, from which we estimate 6o = 218 ± 19 km s~^ for Rq = 8 kpc. 

Our lower limit for the peculiar motion of Sgr A* of about 20 km implies a lower limit for mass of 
Sgr A* of ~ 10^ Mq . This rules out the possibility that Sgr A* is any known multiple star system, such as 
a contact binary containing ^ 10 Mq and radiating near its Eddington luminosity. A mass of more than 
10^ Mq and a luminosity ^ 10^ Lq indicates that Sgr A* is radiating at ^ 0.1% of its Eddington limit. 
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All observations are consistent with Sgr A* being a super-massive black hole. Since the lower limit 
for the mass of Sgr A* is only about 0.1% of the gravitational mass inferred from the stellar motions, 
one cannot claim from our observations alone that even a significant fraction of the dark mass must be 
in Sgr A*. However, alternative models involving "dark" matter distributions are severely restricted by 
observations. 

Future VLBA observations should be able to reduce the uncertainty in the measurement of the motion 
of Sgr A* out of the Galactic plane to ^ 0.2 km s^^ , at which point knowledge of the Solar Motion may 
become the limiting factor. Should the peculiar motion of Sgr A* be less than 0.2 km s~^ , then its mass 
almost certainly exceeds 10^ Mq. Such a large mass tied directly to the radio source, whose size is ;^ 1 AU 
from VLBI observations, would be compelling evidence that Sgr A* is a super-massive black hole. 

We thank D. Backer for providing coordinates for sources J 1745 -283 and J 1748-291 early in our 
project, M. Eubanks for measuring the astrometric position of J1745-283, and V. Dhawan for helping with 
the VLBA setup. 
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Table 1. Residual Position Offsets Relative to Sgr A' 



Source 


Date of 
Observation 


East Offset" 
(mas) 


North Offset" 
(mas) 


J 1745-283 


1995 March 4 


0.73 


1.18 




1996 March 20 


3.60 


5.80 




1996 March 31 


3.75 


6.64 




1997 March 16 


6.90 


11.18 




1997 March 27 


7.10 


11.13 


J1748-291 


1996 March 20 


27.55 


82.07 




1996 March 31 


27.53 


81.98 




1997 March 16 


31.00 


86.92 




1997 March 27 


31.09 


86.60 



'^Uncertainties are 0.1 and 0.4 mas for the East and North 
offsets, respectively, for all epochs except 1995 March 4 which are 
estimated to be 0.5 and 0.8 mas owing to poor phase coherence. 

Note. — Positions are relative to Sgr A*, after removing the 
wO.7 degree differences. The origins are based on our originally 
adopted J2000 positions for Sgr A* (17 45 40.0500, -29 00 28.120), 
J1745 283 (17 45 52.5056, -28 20 26.302), and J1748-291 (17 48 
45.6930, -29 07 39.488). 
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Table 2. Apparent Relative Motions 



Source - Reference 


East Offset 


North Offset 


Offset 


b^^ Offset 




(mas y-^) 


(mas y-^) 


(mas y-^) 


(mas y-^) 


Sgr A* - J1745-283 


-3.33 (0.10) 


-4.94 (0.40) 


-5.90 (0.35) 


+0.20 (0.30) 


J1748-291 - J1745-283 


+0.17 (0.14) 


-0.22 (0.56) 


-0.10 (0.56) 


-0.26 (0.42) 



Note. — Motions are based on data in Table 1, and estimated uncertainties are given in 
parentheses. 



Table 3. Sgr A*'s Motion in Galactic Coordinates" 



Description 




b" 




(km s-i) 


(km s~^) 


Observed Sgr A* motion 


-223 (19) 


8(11) 


Effects of Solar Motion'' removed 


-218 (19) 


15 (11) 


Effects of Galactic Rotation'' removed 


0(15) 


15 (11) 



Motions arc for Sgr A*, based on the J 1745-283 (relative to 
Sgr A*) data in Table 2. Speeds assume i?o = 8.0 ± 0.5 kpc (Reid 
1993). Quoted uncertainties given in parenthesis are 1-cr and include 
measurement uncertainty and an angular-to-linear motion scaling 
error from the uncertainty in i?o • 

''Adopted Solar Motion with respect to a circular orbit is 5.25+0.62 
km in and 7.17+0.38 km in tf' (Dchncn & Binncy 1998). 

Adopted circular rotation of 27.19 + 0.87 km s^^ kpc^^ (Feast 
& Whitelock 1997) removed from our measured angular rotation rate 
of —27.2+ 1.9 km s"'^ kpc""*^ (after correction for the Solar Motion) 
and then multiplied by i?o = 8.0 kpc. Since the angular motion 
difference is small, it is not sensitive to the value adopted for Rq when 
converting to linear motion. The quoted uncertainty is dominated by 
measurement uncertainties for the angular motions, scaled by Rq . 




1995 1996 1997 

Time (years) 

Fig. 1. Position residuals of Sgr A* relative to J1745-283 (circles) and J1748-291 relative 
to J1745-283 (crosses) versus time. Eastward components are shown in the top panel and Northward 
components in the bottom panel. The solid and dashed lines give the variance-weighted best fit components 
of proper motion. The J 1748-291 — J 1745-283 positions have been offset to fit the plot scale for the 
Sgr A*— J1745-283 data. 
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-5 -10 



East Offset (mas) 

Fig. 2. Position residuals of Sgr A* relative to J1745 283 on the plane of the sky. North is to 
the top and East to the left. Each measurement is indicated with an ellipse, approximating the apparent, 
scatter broadened size of Sgr A* at 43 GHz, the date of observation, and 1 — a error bars. The dashed line is 
the variance-weighted best-fit proper motion, and the solid line gives the orientation of the Galactic plane. 
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Fig. 3. Enclosed mass versus radius for various model distributions of dark matter, assuming the 
mass of Sgr A* <C 10® Mq . Models labelled "a" through "d" have decreasing central mass condensations 
(progressively softer gravitational potentials) and approximately bound mass distributions that are 
consistent with stellar proper motion data. Model "a" is a point mass; model "b" through "d" have 
Phimmer density distributions with po of 3.9 x 10^^,2.5 x 10",and 6.0 x 10" Mq pc'^; ro of 0.00002, 
0.001, and 0.01 pc; and a of 3, 4, and 5, respectively. The sloping solid line indicates the upper limit 
for enclosed mass based on the proper motion of Sgr A*. The vertical solid line at 3 x 10"^ pc (6 AU) 
indicates the upper limit in radius, where angular excursions of Sgr A* of < 0.8 mas could be missed 
owing to insufficient astrometric accuracy. The sloping dotted line indicates expected improvement in the 
measurement of the proper motion of Sgr A* within 5 years. 



